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ABSTRACT 
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I/"") | In proto-neutron stars with strong magnetic fields, the cross section for v e {v e ) 

absorption on neutrons (protons) depends on the local magnetic field strength due to 



the quantization of energy levels for the e~ (e + ) produced in the final state. If the 
neutron star possesses an asymmetric magnetic field topology in the sense that the 
magnitude of magnetic field in the north pole is different from that in the south pole, 
^ cj then asymmetric neutrino emission may be generated. We calculate the absorption 

cross sections of v e and v e in strong magnetic fields as a function of the neutrino 
energy. These cross sections exhibit oscillatory behaviors which occur because new 
Landau levels for the e~ (e + ) become accessible as the neutrino energy increases. By 
evaluating the appropriately averaged neutrino opacities, we demonstrate that the 
change in the local neutrino flux due to the modified opacities is rather small. To 
^ . generate appreciable kick velocity (~ 300 km s _1 ) to the newly-formed neutron star, 

the difference in the field strengths at the two opposite poles of the star must be at 
least 10 16 G. We also consider the magnetic field effect on the spectral neutrino energy 
fluxes. The oscillatory features in the absorption opacities give rise to modulations in 
the emergent spectra of v e and v e . 



Subject headings: stars: neutron - pulsars: general - supernovae: general - dense 
matter - magnetic fields - radiation transfer 



1. Introduction 



It has been recognized shortly after the discovery of pulsars that neutron stars have velocities 
much in excess of those of any other normal stellar populations in our Galaxy (Gunn &; Ostriker 
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1970; Minkowski 1970). However, it is only in the last few years that a significant body of evidence 
has come into place to support the view that type II supernovae are nonspherical and neutron 
stars receive large kick velocities at birth. Lyne and Lorimer (1994) analyzed pulsar velocities in 
light of new proper-motion measurements (Harrison et al. 1993) and up-to-date pulsar distance 
scale (Taylor & Cordes 1993), and concluded that pulsars were born with a mean speed of 
~ 450 km s _1 , much larger than previously thought. More recent studies of pulsar velocities have 
adopted more sophisticated satistical methods and included better treatment of selection effects 
and uncertainties (Lorimer et al. 1997; Hansen & Phinney 1997; Cordes & Chernoff 1997). They all 
yielded a result for the pulsar birth velocity in qualitative agreement (although not in quantitative 
details) with that of Lyne & Lorimer (1994). In particular, Cordes and Chernoff (1997) found that 
the three-dimensional space velocities of their sample of 47 pulsars have a bimodal distribution, 
with characteristic speeds of 180 km s -1 and 700 km s _1 (corresponding to 80% and 20% of the 
population). They also estimated that pulsars with velocities greater than 1000 km s _1 may be 
underrepresented owing to selection effects in pulsar surveys [The uncertainty in the high-velocity 
end of the pulsar velocity distribution function has also be emphasized by Hansen & Phinney 
(1997)]. Concrete evidence for the existence of pulsars with velocities of ^ 1000 km s _1 has come 
from the observation of the Guitar Nebula pulsar (B2224+65), which produces a bow shock when 
plowing through the interstellar medium (Cordes, Romani & Lundgren 1993). In addition, studies 
of pulsar-supernova remnant associations have uncovered a number of pulsars having velocities 
much greater than 1000 km s _1 (Frail et al. 1994), although in many cases the associations are 
not completely secure. 

Compelling evidence for supernova asymmetry and pulsar kicks also comes from the detection 
of geodetic precession in the binary pulsar PSR 1913+16 (Cordes et al. 1990; Arzoumanian et 
al. 1996), and the orbital plane precession in the PSR J0045-7319/B star binary (Lai et al. 1995; 
Kaspi et al. 1996) and its fast orbital decay (which indicates retrograde rotation of the B star with 
respect to the orbit; see Lai 1996). These results demonstrate that binary break-up (as originally 
suggested by Gott, Gunn & Ostriker 1970; see Iben & Tutukov 1996) can not be solely responsible 
for the observed pulsar velocities, and that natal kicks are required. In addition, evolutionary 
studies of neutron star binary population imply the existence of pulsar kicks (e.g., Deway &: 
Cordes 1987; Fryer & Kalogera 1997; Fryer et al. 1998; see also Brandt & Podsiadlowski 1995). 
Finally, there are a large number of direct observations of nearby supernovae (e.g., Cropper et 
al. 1988; Trammell et al. 1993; McCray 1993; Utrobin et al. 1995) and supernova remnants (e.g., 
Morse, Winkler & Kirshner 1995; Aschenbach et al. 1995) in radio, optical, and X-ray bands which 
support the notion that supernova explosions are not spherically symmetric. 

The origin of the pulsar velocities is still unknown. Two classes of mechanisms for the natal 
kicks have been suggested^. The first class relies on convective instabilities in the collapsed stellar 



1 These exclude the slow, post-explosion rocket effect due to electromagnetic radiation from off-centered magnetic 
dipole moment (Harrison & Tademaru 1975). 
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core and within the rebounding shock (e.g., Burrows & Fryxell 1992; Burrows et al. 1995; Janka 
& Miiller 1994, 1996; Herant et al. 1994). The asymmetries in the matter and temperature 
distributions associated with the instabilities naturally lead to asymmetric matter ejection and/or 
asymmetric neutrino emission. Numerical simulations indicate that the local, post-collapse 
instabilities are not adequate to account for kick velocities higher than a few hundred km s _1 
(Janka & Miiller 1994). A variant of this class of models therefore relies on the global asymmetric 
perturbations seeded in the presupernova cores (Goldreich et al. 1996; see also Bazan &: Arnett 
1994). Clearly, the magnitude of kick velocity depends on the degree of initial asymmetry in the 
imploding core (Burrows & Hayes 1996). Due to various uncertainties in the presupernova stellar 
models, it is not clear at this point whether sufficiently large initial perturbations can be produced 
in the precollapse core (Lai &: Goldreich 1998). 

In this paper we focus on the second class of models in which the large pulsar kick velocities 
arise from asymmetric neutrino emission induced by strong magnetic fields. Since 99% of the 
neutron star binding energy (a few times 10 53 erg) is released in neutrinos, tapping the neutrino 
energy would appear to be an efficient means to kick the newly-formed neutron star. Magnetic 
fields are naturally invoked to break the spherical symmetry in neutrino emission. But the actual 
mechanism is unclear. We first briefly comment on previous works on this and related subjects. 

1.1. Previous Works on Asymmetric Neutrino Emission Induced by Magnetic Fields 

A number of authors have noted that parity violation in weak interactions may lead to 
asymmetric neutrino emission from proto-neutron stars (Chugai 1984; Dorofeev et al. 1985; 
Vilenkin 1995; Horowitz &; Piekarewicz 1997). However, their studies are largely unsatisfactory for 
a number of reasons: they either failed to identify the most relevant neutrino emission/interaction 
processes or the relevant physical conditions in proto-neutron stars, or stopped at estimating 
the magnetic field effects on neutrino opacities. Chugai (1984) and Vilenkin (1995) (see also 
Bezchastnov & Haensel 1996) considered neutrino-electron scattering and concluded that the effect 
is extremely smallf] (e.g., to obtain Vkick = 300 km s _1 would require a magnetic field of at least 
10 16 G). However, neutrino-electron scattering is less important than neutrino-nucleon scattering 
in determining the characteristics of neutrino transport in proto-neutron stars. Similarly, Dorofeev 
et al. (1985) considered neutrino emission by Urea processes in strong magnetic fields. But as 
shown by Lai & Qian (1998), in the bulk interior of the neutron star, the asymmetry associated 
with neutrino emission is cancelled by the asymmetry associated with neutrino absorption. 

Concerning the parity violation effect in proto-neutron stars, Horowitz & Li (1997) recently 
suggested that the asymmetry in neutrino emission may be enhanced due to multiple scatterings 



2 Note that Chugai's estimate for the electron polarization in the relativistic and degenerate regime (the relevant 
physical regime) is incorrect. This error leads to an overestimate of the effect as compared to Vilenkin's result. 
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of neutrinos by nucleons which are slightly polarized by the magnetic field. They estimated that 
a field strength of a few times 10 12 G is adequate to account for kick velocities of a few hundred 
km s _1 . However, the paper by Horowitz & Li (1997) only discusses the idealized situations of 
scattering media, and ignores various essential physics that is needed in a proper treatment of 
neutrino transport in proto-neutron stars. In particular, it does not consider the effect of neutrino 
absorption which one might suspect to wash out the cumulative effect from multiple scatterings. 
Our own study of the parity violation effect (Lai & Qian 1998) was flawed in the treatment of 
scattering terms in the neutrino transport equation. As shown by Arras & Lai (1998), detailed 
balance requires that there be no cumulative effect from multiple scatterings in the bulk interior of 
the proto-neutron star where local thermodynamic equilibrium applies to a good approximation. 
Enhancement of neutrino emission asymmetry from multiple scatterings obtains only after 
neutrinos thermally decouple from proto-neutron star matter, and therefore is insignificant. 

Bisnovatyi-Kogan (1993) attributed pulsar kicks to asymmetric magnetic field distribution in 
proto-neutron stars. Using the fact that neutron decay rate can be modified by the magnetic field, 
he inferred that neutrino emissions from opposite sides of the neutron star surface are different. 
However, neutron decay is not directly relevant for neutrino emission from a newly- formed neutron 
star. Roulet (1997) (whose paper was posted while this paper was being written) considered the 
relevant neutrino absorption processes. But his study was restricted to calculating the neutrino 
cross sections as functions of neutrino energy (corresponding to §2.1 of this paper), and therefore 
was insufficient for addressing the issue of asymmetric neutrino emission. 

Finally, it is worthwhile to mention a speculative idea on pulsar kicks which relies on 
nonstandard neutrino physics. Kusenko & Segre (1996) suggested that asymmetric v T emission 
could result from the Mikheyev-Smirnov-Wolfenstein flavor transformation between v T and v e 
inside a magnetized proto-neutron star because a magnetic field changes the resonance condition 
for the flavor transformation. With a v T mass of ~ 100 eV, they claimed that magnetic fields 
of ~ 3 x 10 14 G can give the pulsar a kick velocity of a few hundred km s^ 1 . However, their 
treatment of neutrino transport was oversimplified: e.g., they ignored that neutrinos of different 
energies have different resonance surfaces. Furthermore, a simple geometric effect and realistic 
proto-neutron star conditions easily reduce their estimated kick velocity by an order of magnitude. 
Most likely, a magnetic field strength of ~ 10 16 G is needed to produce the observed average pulsar 
kick velocity via their mechanism (see Qian 1997). Similarly strong magnetic fields are required 
in variants of the Kusenko & Segre mechanism, e.g., those considered by Akhmedov et al. (1997) 
(which relies on both the neutrino mass and the neutrino magnetic moment to facilitate the flavor 
transformation) . 

1.2. Plan of This Paper 

As should be clear from the brief review in §1.1, in spite of quite a number of suggestions, 
there is at present no consensus on the magnitudes of the magnetic field induced asymmetry in 
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neutrino emission from proto-neutron stars and the resulting kick velocities. In this paper we study 
the effect of asymmetric magnetic field topology on pulsar kicks. Since the energy levels of e~ and 
e + in a magnetic field are quantized, the v e and v e absorption opacities near the neutrinosphere 
and the neutrino-matter decoupling region depend on the local magnetic field strength. If the 
magnitude of magnetic field in the north pole is different from that in the south pole (the field 
does not need to be ordered), then asymmetric neutrino flux may be generated. Here we ignore 
the effects of magnetic fields on the proto-neutron star structure through the equation of state - 
these are secondary effects as far as neutrino transport is concerned. 

The organization of this paper is as follows. In §2 we calculate the v e and v e absorption cross 
sections as functions of neutrino energy in strong magnetic fields. Section 3 summarizes the basic 
features of neutrino transport near the neutron star surface. In §4 we derive and evaluate the 
"Rosseland mean" neutrino opacities which are directly related to the local neutrino flux. We 
then consider in §5 the change in neutrino flux due to the modified opacities, and estimate the 
kick velocity resulting from an asymmetric magnetic field topology. In §6 we briefly consider how 
the emergent neutrino spectra may be modified by the strong magnetic field. In §7 we present our 
conclusion. 

Unless noted otherwise, we shall use units in which h, c, and the Boltzmann constant k are 
unity. 



2. v e and v e Absorption Cross Sections in Strong Magnetic Fields 

In this section we consider how magnetic fields modify the (total) neutrino absorption 
{y e + n — > p + e~ and U e + p — > n + e + ) cross sections for specific neutrino energies. We present 
our calculation of the v e absorption cross section in detail. The 9 e absorption cross section can be 
obtained in a similar manner. 

To begin with, it is useful to review various energy scales involved in the problem. Near the 
neutrinosphere and the neutrino-matter decoupling region, the matter density is typically 10 11 -10 12 
gem -3 , and the temperature is about several MeV. The electrons are extremely relativistic and 
highly degenerate, with Fermi energy Ep, e — 51.5 (Yep^) 1 ^ 3 MeV, where Y e ~ 0.1 is the electron 
fraction (the number of electrons per nucleon) and p±2 = p/(W 12 gem -3 ). The neutrons and 
protons have Fermi energies Ep, n = 1-4 [(1 — Y e )pi 2 ] 2 ^ 3 MeV and Ep, p = 1.4 (Y" e pi 2 ) 2 / 3 MeV, 
respectively, and thus they are essentially nondegenerate. 
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2.1. The Cross Sections 

In strong magnetic fields the transverse motion of the electron is quantized into Landau levels. 
The electron energy dispersion relation is 



E e = 



p\ + m\{\ + 2nb)} 1/2 , n = 0, 1, 2, ••• (1) 



where p z is the longitudinal momentum (along the magnetic field), and n is the Landau level 
index. The dimensionless field strength b is defined as 

b = —, with B c = = 4.414 x 10 13 G. (2) 

B c he 

[B c is the critical field strength defined by equating the cyclotron energy heB / (m e c) to m e c 2 .] The 
quantization effect on the proton is extremely small due to the large proton mass, and is neglected 
throughout this paper. 

Since the cyclotron radius (characterizing the size of the Landau wavefunction) 
\hc/(eB)] 1 / 2 = [h/(m e c)]b~ l l 2 = 386 b~ x l 2 fm is much greater than the range of the weak 
interaction [h/(mwc) ~ 2.5 x 10~ 3 fm, where mw — 80 GeV is the mass of W boson], one can 
use the V — A theory with contact interaction. Locally (at the point of interaction) the Landau 
wavefunction resembles a plane wave. Thus we expect that the spin-averaged matrix element is 
the same as in the case without magnetic fields. This expectation has been borne out by detailed 
calculations (e.g., Fassio-Canuto 1969; Matese & O'Connell 1969). Therefore we only need to 
consider the effect of magnetic field on the electron phase space. 

Since the neutrino energy (a few to tens of MeV) of interest is much less than the nucleon 
mass, we first consider the case where the nucleon mass is taken to be infinity, an approximation 
adopted in the standard zero- field calculations (e.g., Tubbs k, Schramm 1975). The cross section 
for v e + n — > p + e~ can be written as a sum over the (spatial and spin) states of the electron (per 
unit volume) 

e 

where we have defined 

A = ttG 2 f cos 2 c (cy + 3<£), (4) 

and the other symbols have their usual meanings: Gp = (293 GeV) -2 is the universal Fermi 
constant, cos 2 6c = 0.95 refers to the Cabbibo angle, cy = 1 and ca = 1-26 are the coupling 
constants, Q = m n — m p = 1.293 MeV is the neutron-proton mass difference, E u is the v e energy, 
and (1 — f e ) is the Pauli blocking factor. The electron occupation number f e is 

/e= exp [{E e -» e )/T\ + V (5) 
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where E e is the electron energy, \i e is the electron chemical potential, and T is the matter 
temperature. For zero field, the sum becomes 2 J d 3 p/(27r) 3 , and we have 



4 a =o(^) = + Q) ( E « + Q) 2 - < "' (i - A), ( 6 ) 



.4 



7T 



1/2 



with E e = E u + Q. For B ^ 0, knowing that the degeneracy of the Landau level (per unit area) is 
g n eB/hc = g n m 2 b / (2ir) , with g n = 1 for n = and g n = 2 for n > 0, we have 



4 ab8) (^) = A f°° £ 5 n^(l - f e ) 5{E U + Q-E e ) = A^4 f> n ^(l - f e ), (7) 
J-oo 2tt ^ 2vr 2vH ^ p z , n 

where £ e = E v + Q, and p z , n = [E 2 - m 2 e (l + 2n6)] 1/2 . The upper limit n max for the sum is the 
maximum value of n for which p z n is meaningfully defined, i.e., 



^max — lilt 



{Eu + Qf-m 2 
2m 2 b 



(8) 



where Int[x] stands for the integral part of x. Henceforth, n max is similarly defined, and we shall 
not explicitly write out the lower and upper limits for the sum. 

Figure 1 depicts ^^^{Ey) for b = 100 and a^^{E v ) as functions of E v for typical conditions 
(Y e = 0.1, /i e = 20 MeV, and T = 3 MeV) in a proto-neutron star. Note that ideally, one should 
compare a^° B \E v ) and (t^^q(E u ) at the same matter density. For /j, e ~ Ep, e 3> T (so that 
n e - 3> n e +), the matter density is related to \i e via 

= ^iJ^ f><We, (for 5 = 0), (9) 

P = T^2^ 



/•oo 

/ d Pz Y,9nfe, (for B^0). (10) 



In practice, we have found that for a given /x e , the dependence of /? on the field strength is 
rather weak for the range of parameters of interest in this paper. Thus in all of our calculations 
here and below, we compare the cross sections at a given /x e . A prominent feature of a^ bs \E v ) 
shown in Fig. 1 is the presence of "resonances," which occur when p z<n = (cf. Eq. J7j), i.e., at 
the neutrino energies for which the final state electron lies in a "stationary" Landau level (with 
zero longitudinal momentum). These resonances reflect the divergence of the density of states of 
electrons at these stationary Landau levels. 

The cross section for v e + p — > n + e + can be similarly calculated. One only needs to replace 
Q with —Q and \x e with — ju e in the expressions for v e absorption to obtain the desired expressions 
for v e absorption. The numerical results for v e absorption cross section are shown in Fig. 2 for the 
same physical parameters as used in Fig. 1. 
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2.2. Thermal Averaging 

The divergent cross sections at the "resonances" shown in Figs. 1 and 2 are clearly unphysical. 
They are smoothed out by the thermal motion of the neutron. To the leading order, the width of 
the resonance is (fcT/mn) 1 / 2 ^, corresponding to the Doppler shift of the neutrino energy in the 
neutron's rest frame. The recoil energy, of order E^/m n , is much less effective in smoothing out 
the resonances, and will be neglected. With this approximation, the smoothed cross section can 
be written as 

(4 abs) (^)} ~ — )2 jvldv n f n jdco S 6 n a { t\K), (11) 

where E' v = E u (l — v n cos 6 n ) and /„ = (27rm n /T) 3 / 2 exp( —m n v^/2T) is the velocity distribution 
function for the neutron. If we replace the factor (1 — f e ) [which depends on (E' v + Q)] in 
(jg hs \El) with the average value (1 — f e ), which can be approximately taken as (1 — f e ) evaluated 
at E e = E v + Q, then the angular integral in Eq. ([ll]) can be carried out analytically. For B = 0, 
we have 



(4 a = ) (^)) 

which is reduced to 



A 



f2vr) 2 



(1 " fe) J VndVn / n ^-(^ - m 2 ) 3 / 2 



E e =E„(l+v n )+Q 



E e =E v (l-v n )+C 



(12) 



(a^(E u ))c^(E v + Qf(l-f e 
if we neglect the electron mass m e . For B ^ 0, we find0 



1 + 



2 rp 

E V + QJ m,„ 



(13) 



(ef s) {E u )) 



(1 - f e ) IVndVnfn- 



(2VT) 



! 2vr 2 K 



9n [E 2 e - m 2 (l + 2nb) 



1/2 



E B =E v (l+v n )+Q 



E e =E v (l-v n )- 



(14) 



We could have approximately included the neutron recoil effect by replacing Q with 
Q — (E v + Q) 2 /(2m n ) in Eqs. (fi^2|)-(|r3). But we choose to leave out this effect since as 
mentioned previously, it is ineffective in smoothing out the resonances in s (E u ). We have 
numerically evaluated Eqs. (|l^) and (|l4]), and the results are depicted in Fig. 1 together with 
the cross sections before averaging. Similar results for the v e absorption cross section are shown 
in Fig. 2. As expected, we find that for B = 0, the averaged cross section (&b=o(Ev)) 1S almost 
identical to a^^(E v ) (see Eq. [13]). However, for B / 0, the thermal averaging is essential in 



obtaining the smoothed physical absorption cross sections. 



3 To obtain this result, one has to interchange the sum and the angular integral. This interchange is valid because 



\l/2 



has an integrable singularity at x - 
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2.3. Limiting Cases 



The oscillatory behaviors of (a^g 38 ' (E u )) in Figs. 1 and 2 can be demonstrated analytically. 
Using the Poisson summation formula (see e.g., Ziman 1979): 



OO -i OO /*oo 

E/fr+gH £ ("W f(x)e 2msx dx, 

n=0 ~ s=—oo " 



(15) 



we can rewrite the sum in Eq. ( |14|) as 

m 2 e bY,9n [E 2 e -m 2 e {l + 2nb) 



1/2 



1 



s=l 



exp 



2iris 



Ei 



2m 2 b 



(16) 



The first term on the right-hand side of Eq. (|i~6| ) corresponds to the zero-field result, and the 
second term (involving the sum of many oscillatory terms) corresponds to the correction due to 
the quantization effect. The sum is clearly convergent, and thus for small B, the averaged cross 
section (a^ bs \E u )) is reduced to the zero-field value (a^^j (E u )} . 

For sufficiently large B at a given E u , or, equivalently, for sufficiently small E u at a given B, 
only the ground state Landau level is filled by the electron, i.e., n max = (see Eq. ||). In this 
limit, {&g (E v )) > {E u )) /4, and only the first term in Eq. (|7|) needs to be retained. We 

then have 

mlb (E u + Q) _ 

211/2 ^ !e > 



(4 abs) 



(E V ))~A. 



(17) 



2vr 2 [{E u + QY-m 

for v e absorption. A similar expression for v e absorption is obtained by replacing Q with — Q and 
/j, e with — [i e in Eq. ([17]). 



2.4. Discussion 

The physical conditions which Figs. 1 and 2 are based on (T = 3 MeV, Y e = 0.1, and 
He = 20 MeV, corresponding to a matter density of p = 7.2 x 10 11 g cm -3 ) are typical of the 
neutrinosphere of the proto-neutron star. In Fig. 3 we plot the v e and v e absorption cross sections 
for a number of slightly different physical parameters, specified by b, \x e and T (all with Y e = 0.1). 
The cross section for elastic neutral-current scattering v + ./V — > v + N, a^ sc \E v ) ~ G 2 pE 2 j-K 
(e.g., Tubbs & Schramm 1975), is also shown for comparison. Absorption dominates the transport 
of f e for high energies, with the ratio of absorption to scattering cross section per nucleoli 
(1 — y e )a"( abs ) / 'o"( sc ) approaching 5 for E u 3> /U e . For smaller energies {E v <> /i e ), absorption is 
severely suppressed by the Pauli blocking of the final electron states, and scattering becomes 
more important. In the case of v e , Pauli blocking is negligible, and scattering almost always 
dominates the transport, with the ratio of scattering to absorption cross section per nucleon 
tf( sc V[>> (abs) ] oi 2(0.1/y e ) for Eu ^ 10 MeV. 
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As we can see from Figs. 1-3, for a given set of physical parameters of the medium (T, Y e , fj, e 
or p), the magnetic field effect on the absorption opacity is more prominent for lower-energy 
neutrinos than for higher-energy ones. This is because for a given field strength, as the neutrino 
energy increases, more and more Landau levels become available for the electron in the final state 
(See Eq. ||). Consequently, the electron phase space and hence the cross section should approach 
the zero-field values according to the correspondence principle. We found that typically when the 
number of Landau levels filled by the electron (n max ) is greater than 2-3, the magnetic field effect 
on the absorption opacity becomes negligible. The "critical" neutrino energy is approximately 
given by 

4 crit) - m e (26n max + 1)^ * 6 (^p) ^ (^) ^ MeV. (18) 

For E v <; ii£ c , one can approximate (a^ (E„)) by the zero-field value {ct^^q {E v )) with good 
accuracy (within a few percent). Alternatively, we can define a threshold magnetic field, B t ^, 
above which the v e and v e absorption opacities are modified (by more than a few percent) by the 
magnetic field: 

B -S B - 3 * iol5 (i^) 2 Gi) G - (i9) 

where (E v ) ~ 10 MeV is the average neutrino energy. 

In the deeper layers of the neutron star, the typical neutrino energy is higher. Therefore 
we expect the magnetic field effect on the absorption opacities to be smaller there. In addition, 
for f e , absorption becomes unimportant as compared to scattering in the high density region 
(although not necessarily in the deepest interior, where v e can be degenerate). Thus, in order to 
assess the importance of the magnetic field effect on the neutrino emission, we only need to focus 
on the surface region, the region between the neutrinosphere and the neutrino-matter decoupling 
layer. We note that because the absorption opacities in the magnetic field for different neutrino 
energies oscillate around the the zero-field values, without a detailed calculation of the mean 
opacity (appropriately averaged over the neutrino energy spectrum) one cannot even answer the 
qualitative question such as whether the magnetic field increases or decreases the local neutrino 
energy flux. 

In §§3-6, we study how the modified absorption cross sections affect the emergent v e and v e 
energy fluxes and spectra from the proto-neutron star. 



3. Neutrino Transport Near the Stellar Surface 

As discussed in §2.4, the quantization effect of magnetic fields on neutrino absorption is 
important in the region only near the neutron star surface (between the neutrinosphere and the 
neutrino-matter decoupling layer). In this section we consider neutrino transport in this relatively 
low density region. 
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Because u e and v e with typical energies of 10 MeV have comparable absorption and scattering 
opacities (see Fig. 3), the emergent v e and v e energy fluxes and spectra are determined essentially 
in the same region where these two neutrino species just begin to thermally decouple from 
the proto-neutron star matter. The v e and v e occupation numbers in this region are then 
approximately given by 

1 

exp[(E v /T) -rj u \ + l 

and 

fu e \E u ) = 71— T=7 j : — (21) 

exp[(E u /T) + r] v \ + 1 

respectively, where rj v = p u /T, and \i v is the v e chemical potential. From these occupation 
numbers, we obtain 

^ e -np e = ^fl + 4V 3 , (22) 



and 



6 \ 7T 



U +U- -!!L!(i + M + lMW 4 (23) 



where for example, n Ue and U Ue are the number and energy densities, respectively, for the v e . We 
can estimate the magnitude of r) v as follows. The net deleptonization rate of the proto-neutron 
star is 

• Y eA M/m N 2 
- N e ~ — '— ~ c(n Ue - upJttR , (24) 

Miff 

where Y e ^ is the initial electron fraction of the proto-neutron star, M is the proto-neutron star 
mass, rriN is the nucleon mass, t^is is the neutrino diffusion timescale, and R is the proto-neutron 
star radius. The sum of the v e and v e luminosities is 

1GM 2 /R % , 

3 Cdiff 

From Eqs. (p2|)-(p5|), we have 

217T 2 1 + (30/7^)^ + (15/7vr 4 )^ *^Tg / F e ,, \ / T \ 

* ~ - 20 1 + (r^) GM^ ~ °- 1 lo36 J loWj « X ' (26) 

where the numerical value is obtained for M ~ IAMq and R ~ 10 km. Therefore, to a good 
approximation, the and P e in the decoupling region have Fermi-Dirac energy distributions with 
zero chemical potentials. 

Let Ie = -^E( r )^) be the specific intensity of neutrinos of a given species, where r is the 
spatial position, and ft specifies the direction of propagation. The general transport equation 
takes the form 

n ■ vi E = P ^ hs) *i { E FD) + p4 c) ^ J dn'i' E - P [4 bs) * + 4 C) ] Ie, (27) 
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where I' E = Ie( t , « )' an( ^ the subscript U E" indicates that the relevant physical quantities depend 
on the neutrino energy E. The opacities K E bs ^* and related to the corresponding cross 

sections as, e.g., k^^* = (1 — Y e )<jg bs ^* /mjv and = 0£ / m N for z^ e . Note that we have 
included the effect of stimulated absorption for neutrinos (e.g., Imshennik & Nadezhin 1973) by 
introducing the corrected absorption opacity £ j [l + exp(-£/T)] (T is the matter 

temperature). Hereafter, we shall suppress the superscript "*" in k^^*, and the notation K^ bs ^ 
should be understood as having included the correction due to stimulated absorption. The quantity 
I^ D ^ in Eq. ( f27|) is the equilibrium neutrino intensity, given by the Fermi-Dirac distribution 

{FD) E 3 1 

E c 2 /i 3 exp(£/T) + 1' 1 ' 



Note that the time derivative term (OIe /dt) on the left-hand-side of Eq. (|27| ) has been neglected. 
This corresponds to an instantaneous redistribution of emission and absorption sources (i.e., an 
instantaneous redistribution of matter temperature) — The timescale for the redistribution is of 
order the mean free path divided by c, much smaller than the neutrino cooling time of the star. 

The zeroth order moment of the transport equation can be written as 

V.F £ = p4 bs) c(4 FD) -U E ), (29) 

where Fe = JlE&dtl is the (spectral) neutrino energy flux, Ue = JlEdQ/c is the energy density, 
and Ug' = (Att/c)!^ is the corresponding equilibrium energy density. With the Eddington 



closure relation, the first moment of Eq. (27) becomes 

„ c 



m VU E , (30) 



where = K^ hs ^ + Kg C \ Equations (|29D and (p0) completely specify the neutrino radiation field. 



4. Rosseland Mean Opacities 

Here we are concerned with the total neutrino energy flux emitted from a local surface region 
of the neutron star. The magnetic field strength in this local region is assumed to be constant. 
Since the scale heights for physical quantities such as the matter density and temperature are 
much smaller than the local radii, we adopt the plane-parallel geometry. Integrating Eq. (|30| ) 
over the neutrino energy, we find that the total neutrino energy flux (of a specific species) 
F = fFgdE = Fr is given by 

3 p kr dr 

where U = J UEdE is the total neutrino energy density, and kr is the "Rosseland mean opacity" 
defined as 

1 fdE(l/K$)(dU E /dr) 



kr JdE (dU E /dr) 



(32) 
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Since neutrino emission and absorption rates are extremely fast even in the regions near the 
proto-neutron star surface, a radiative equilibrium^] approximately holds, i.e., 

p K g bs ) (U E FD) -U E )dE~0. (33) 



^From Eq. (29), we have dF/dr ~ 0, or F is approximately constant. 



If we introduce the Rosseland mean optical depth measured from the stellar surface via 
dTR = —pKRdr, with F constant, Eq. (31) can be integrated to yield 



cU = 3F (t r + |) , (34) 

where we have used the boundary condition F = cU/2 at the surface (tr = 0). Equation 
(B3) resembles the standard Eddington profile. Note that outside the decoupling region the 
neutrino energy density U may not be equal to U^ FD "> = (77r 4 /240)T 4 specified by the matter 
temperature T. Similarly, Ue may differ substantially from U^ D \t) (see §6). However, we can 
still approximate Ue by the Fermi-Dirac function at the neutrino temperature T u , which is equal 
to the matter temperature at the decoupling layer. Equation (p%) can then be rewritten as 

k r l^Tl J ^ K W [eME/Tv) + If " W 

With the cross sections given in §2, we can calculate the Rosseland mean opacity k r (B) in 
magnetic fields using Eq. (35). Note that in Eq. ( |35[ ) Kg depends on the local matter temperature 



T. Only in the regions not far from the decoupling layer, do we have T u = T. 

Figure 4 shows the ratio of the Rosseland mean opacity k r (B) of v e to the zero-field value 
k r (0) as a function of the magnetic field strength for typical conditions near the proto-neutron 
star surface. We choose these conditions to be T = T u = 3, 5 MeV, and fj, e = 10, 15, 20 MeV, 
all with Y e = 0.1. Figure 5 shows similar results for v e with similar physical parameters 
(T = T u = 3, 5, 7 MeV; Note that the results are insensitive to fi e ). We see that deviation of 
kr{B) from the zero-field value by more than 5% generally requires field strength b = B/B c J> 100. 

The behavior of the ratio Kr(B) J kr(0) as a function of b is certainly not obvious a priori. 
Since the Rosseland mean integral is dominated by neutrinos with the smallest opacities, and 
since Kg scales linearly with B in the regime where only the ground Landau level is filled (see 
Eq. fi~7f ), we expect that for sufficiently large B, the Rosseland mean opacity increases linearly 
with increasing B. This is indeed what is indicated in Figs. 4 and 5. For "intermediate" field 



4 In principle, the overall radiative equilibrium involves matter and neutrinos of all species. However, because the 
coupling between matter and v T , and v T is much weaker, the overall radiative equilibrium in the region near 

the proto-neutron star surface mainly involves matter and v e and v e . Since the coupling between matter and v e 
is comparable to that between matter and i> e , we further assume that an approximate radiative equilibrium holds 
separately for !/ e and v R . 



-14- 



strength, the Kg s curve for v e (see Fig. 1) has a significant portion which is below the zero- field 
value at low neutrino energies. For such intermediate fields, we would expect kr{B) to be less 
than the zero-field value. This explains the non-monotonic behavior of kr(B) / kr(0) as a function 
of B shown in Fig. 4. For u e , the non-monotonic behavior is much less pronounced (see the curve 
for T = 7 MeV in Fig. 5) since the energy range in which Kg hs ^ < is smaller. 



5. Differential Neutrino Fluxes and Neutron Star Kick ^From Asymmetric 

Magnetic Field Topology 

Imagine a proto- neutron star with an asymmetric magnetic field topology, e.g., the north 
pole and the south pole have different magnitudes of magnetic field (say, by a factor of a few). 
What is the difference in the neutrino fluxes from these two regions due to the different absorption 
opacities of u e and u e , and what is the resulting kick velocity of the neutron star? We recall that 
the fractional asymmetry a in the total neutrino emission (of all species) required to generate a 
kick velocity V kick is a = 0.028 [VW/(10 3 km s" 1 )] [E tot /(3 x 10 53 erg)], where E tot is the total 
neutrino energy radiated from the neutron star (the neutron star mass has been set to 1.4 Mq). 

To quantitatively understand how the magnetically-modified opacities change the local 
neutrino flux, we will have to evaluate these opacities at different depths (corresponding to 
different physical conditions) below the neutron star surface, and to perform complete cooling 
calculations which take into account these opacities. However, Figs. 4 and 5 clearly indicate that, 
even near the stellar surface, where one expects the modification to the opacities to be most 
prominent, and even with the most favorable conditions for the magnetic field to be effective in 
changing the opacities (i.e., low density and low temperature), a 5% change in the Rosseland mean 
opacity (for v e or u e ) requires a magnetic field with b <; 100, or B J> 4 x 10 15 G. In the following, 
we shall content ourselves with merely estimating the change in the local neutrino flux from the 
modified opacities. 

We first consider the case in which the matter temperature throughout the neutron star is 
unaffected by the magnetic field. One might imagine this as a possibility if convection dominates 
energy transfer between different parts of the star. Under the neutrino-matter decoupling layer the 
specific neutrino energy density Ue is the same as the Fermi-Dirac function U^ D \t) specified 
by the matter temperature. The emergent neutrino flux is then inversely proportional to the 



Rosseland mean optical depth tr from the decoupling layer to the stellar surface (see Eq. [34]). 
Since tr for v e or v e at the decoupling layer is of order a few, the flux F is inversely proportional 
to kr evaluated near the decoupling region. If we use the curve for T = 3 MeV and fi e = 10 MeV 
in Fig. 4 (which has the most pronounced magnetic correction), we find that the fractional change 
in the v e flux is 5F Ue /F Ue ~ 0.1 at b = 100 (If we used the curve for \i e = 15 or 20 MeV, which 
more closely represents the physical condition at the decoupling region, the modification 5F Ue /F Vf , 
would be smaller by a factor of a few or more) . Similarly, using the curve for T = 3 MeV in Fig. 5 
to maximize the magnetic correction, we find 8Fp e /Fp e ~ —0.02 at b = 100. Since the fluxes of 
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other types of neutrinos are unchanged^], and since different neutrino species more or less carry 
away the same amount of energy, the change in the total neutrino flux is only 5F to t /-Plot ~ 0.01 
at b = 100. If we now imagine that the north pole of the neutron star has b ~ 100 while the 
south pole has a field strength a factor of a few smaller, then, taking into account the geometric 
reduction (~ 1/3), the flux asymmetry is about a ~ 0.003, just enough to give a kick velocity of 
~ 100 km s _1 to the neutron star. With kr evaluated at more realistic conditions, we expect the 
resulting a to be even smaller. 

The assumption that the matter temperature is unaffected by the magnetic field (and 
therefore spherically symmetric) is unlikely to be valid. In reality, when a region of the star (say 
the north pole) emits more flux, it also cools faster. With reduced temperature, this region will 
tend to emit less flux. In the extreme case, one would find F ~ cUq/tr, with Uq the neutrino 
energy density in the stellar core, and tr the Rosseland mean optical depth from the core to the 
surface. Since in the bulk interior of the star, modification to kr due to the magnetic field is 
negligible, we would conclude that the change in the neutrino flux is much smaller than what is 
estimated in the last paragraph. 

To summarize, the modifications to the absorption opacities of v e and v e due to quantized 
Landau levels give rise to rather small change in neutrino fluxes from the neutron star. To generate 
appreciable kick velocities (~ 300 km s^ 1 ) from this effect requires the difference of field strengths 
at the two opposite poles to be at least 10 16 G (and possibly much larger). 



6. Magnetic Field Effect on v e and v e Energy Spectra 

We now consider how the oscillatory behaviors (see Figs. 1-3) of the absorption opacities of 
Vf, and v e may manifest themselves in the radiated neutrino energy spectra. For an absorption 
dominated medium, neutrinos with energy E decouple from the matter and start free-streaming 
at a column depth (in g cm~ 2 ) y = Jpdr ~ 1/Kg bs \ Thus the neutrino energy spectrum is given 
by Fe ~ cUj^ with evaluated at the depth y ~ I/k^^. In the presence of significant 

scattering, decoupling occurs at depth y ~ l/«^ with the effective opacity K e g - • ■ " 



whereas the neutrinosphere from which neutrinos start free-streaming is located at y ~ I/k^). 
Near the decoupling region we have dUE/dr ~ —pK ( gU E FD \ and thus from Eq. (30) the neutrino 
spectral energy flux is given by 

eff 

FE-^cur\y-V^)- (36) 



5 The decoupling layers of tV( T ) and v„( T ) depend on the v-e scattering cross sections, which can be affected by 
the magnetic field. However, the layers are located much deeper than those of v a and v e (since v-e scattering opacity 
is small), the modification to the v-e scattering opacity due to electron Landau levels is therefore smaller than the 
modification to the absorption opacities of v^ and v E considered in this paper. 
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The above qualitative argument can be made more rigorous by directly solving Eqs. (|29|) and 
]). Consider a plane-parallel atmosphere and define the (energy-dependent) optical depth via 
cLte = —pKgdr = k e dy. Equations (^) and ([?(]) can be combined to give 

,2 (abs) 

U E = (5l[U E -U ( i D \ with pl = 3lE. (37) 

ar E K E ' 

With a constant (3e, the above equation can be solved exactly. For the boundary condition 
Fe = {c/3)(cIUe /cIte) = cUe/2 at te = 0, the neutrino energy flux at the surface is given by 

F E = 1 + ^ /3 fdTEPEClI^iTE) eM-PETE) ^^3 ^^^ * (38) 

Figure 6 shows the spectral energy fluxes of v e calculated using Eq. (|3^) for different field 
strengths (b = 100 and 300). The matter temperature profile is obtained from Eq. ( p4| ) assuming 
U oc T 4 . For a total energy flux of F = (c/A)U(T e g) (where T e g is the effective temperature), 
we have T = T c q[(3 / 4)tr + 1/2] 1 / 4 . We choose T e g- = 3 MeV in these examples. For simplicity, 
we have neglected the variations of the opacities as a function of te- Thus the temperature 
at decoupling (t e ~ see Eq. @) is given by T ~ T cS [(3/4)k r / K { f ] + 1/2} 1 / 4 , where 

Kg ff ' = Pe^e = \[^e^~^e an d we take the relevant opacities to have the corresponding values 
for Y e = 0.1, /j, e = 20 MeV, and T = 3 MeV. This simplification should not introduce qualitative 
changes in the spectrum. As seen from Fig. 6, the spectrum is pinched at the high energies 
because of the -E 2 -dependence of the opacities. The magnetic field introduces modulations in the 
spectrum. For lower magnetic field strengths, the modulations occur at the lower energy band. As 
the magnetic field gets stronger, the modulations shift to the higher energy part of the spectrum. 

The distortions and modulations of the v e and v e energy spectra by strong magnetic fields 
may offer an interesting possibility to limit the magnetic field strength in proto-neutron stars 
through the detection of neutrinos from a Galactic supernova. We note, however, that the u e 
spectra shown in Fig. 6 are based on very approximate calculations. They serve only as illustrative 
examples of the magnetic field effect on the neutrino spectra. More detailed calculations are 
needed to obtain realistic spectra which include the magnetic field effect. 



7. Conclusion 

In this paper, we have studied the issue of whether asymmetric magnetic field topology 
in a proto-neutron star can induce asymmetric neutrino emission from the star through the 
modifications of the neutrino absorption opacities by the magnetic field. These modifications arise 
from the quantized Landau levels of electrons and positrons produced in strong magnetic fields. 
By calculating the appropriate mean neutrino opacities in the magnetic field, we demonstrate that 
this mechanism is rather inefficient in generating a kick to the neutron star: To obtain appreciable 
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kick velocities (~ 300 km s 1 ), the difference in the field strengths at the two opposite poles of the 
star must be at least 10 16 G. 

This work was started while D.L. held the Richard C. Tolman postdoctoral Fellowship 
at Caltech. Additional support was provided by NASA grant NAG 5-2756 and NSF grant 
AST-9417371, as well as a start-up fund for new faculty at Cornell University. Y.Q. is supported 
by the David W. Morrisroe postdoctoral Fellowship at Caltech. 
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Energy of v e (MeV) 



Fig. 1. — Absorption cross section (per nucleon) of u e , (1 — Y e )a^ (E v ), in units of A x MeV 2 = 
ttG 2 f cos 2 8c{cy + 3c\) x MeV 2 ~ 0.91 x 10~ 42 cm 2 , as a function of the neutrino energy. The 
dashed line is the zero-field result (The effect of the nucleon thermal motion is included but is hard 
to discern). The heavy and the light solid lines are the results for b = 100, with and without the 
effect of the nucleon thermal motion, respectively. The dotted line is the elastic scattering cross 
section, a^ sc \E v ), in the same units. All results are for T = 3 MeV, Y e = 0.1, and fj, e = 20 MeV 
(corresponding to p = 7.2 x 10 11 gcm~ 3 ). 
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Fig. 2. — Absorption cross section (per nucleon) of u e , Y e a^ s \E V ). The physical parameters, 
units, and labels are the same as in Fig. 1. 
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Fig. 3. — Absorption cross sections (per nuclcon) of v e (left panel) and v e (right panel), in units 
of A x MeV 2 = ir G 2 F cos 2 9c (cy + 3c 2 A ) x MeV 2 ~ 0.91 x 10~ 42 cm 2 , as functions of the neutrino 
energy for different physical parameters. All results assume Y e = 0.1. The solid, the short-dashed, 
the long-dashed, and the dot-dashed curves are for (b, /x e /MeV, T/MeV) = (100,20,3), (100,10,3), 
(200,20,3), and (100,20,5), respectively. [Note that for v ei the curve for (6, /i e /MeV, T/MeV) = 
(100, 10, 3) is almost identical to the curve for (6, ^ e /MeV, T/MeV) = (100, 20, 3), and both curves 
are close to the curve for (6, /x e /MeV, T/MeV) = (100, 20, 5).] The light dotted lines are the elastic 
scattering cross sections in the same units. 
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Fig. 4. — The ratio of the Rosseland mean opacity kr(B) of u e to its zero-field value kr(0) as a 
function of the magnetic field strength b = B/B c (with B c = 4.414 x 10 13 G). All results assume 
Y e = 0.1. The heavier curves correspond to T = T v = 3 MeV, and the lighter cuvres correspond to 
T = T v = 5 MeV. The solid, the short-dashed, and the long-dashed lines are for fi e = 10, 15, and 
20 MeV, respectively. 
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Fig. 5. — The ratio of the Rosseland mean opacity kr(B) of v e to its zero-field value kr(0) as a 
function of the magnetic field strength b = B/B c (with B c = 4.414 x 10 13 G). All results assume 
Y e = 0.1 and fi e = 20 MeV (but the results are rather insensitive to the value of fi e ). The solid, 
the short-dashed, and the long-dashed line are for T = T u = 3, 5, and 7 MeV, respectively. Note 
that the results mainly depend on T u , and are rather insensitive to T. 
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Fig. 6. — Illustrative examples of the v e spectral energy fluxes (in arbitrary units) modified by the 
magnetic fields. The light solid, the short-dashed, and the long-dashed line are for 6 = 0, 100, and 
300, respectively. The light dotted line is the Fermi-Dirac spectrum with T = T e g- = 3 MeV. 



